Airy functiont 

(exact WKB results for potentials of odd degree) 



Andre Voros 

CEA-Saclay, Service de Physique Theorique 
F-91191 GIF-SUR-YVETTE CEDEX (France) 
(e-mail: voros@spht.saclay.cea.fr) 



An exact WKB treatment of 1-d homogeneous Schrodinger operators (with the 
confining potentials q N , N even) is extended to odd degrees N. The resulting formalism 
is first illustrated theoretically and numerically upon the spectrum of the cubic oscillator 
(potential |(/| 3 ). Concerning the linear potential (N = 1), the theory exhibits a duality 
in which the Airy functions Ai, Ai' become paired with the spectral determinants of 
the quartic oscillator (N = 4). Classic identities for the Airy function, as well as some 
less familiar ones, appear in this new perspective as special cases in a general setting. 

A number of quantum spectral properties have now become established (by means of 
WKB theory, asymptotic or exact) for homogeneous q N potentials on the real line having a 
discrete spectrum (confining case, i.e., N even). After recalling the needed results (Sec. 1), 
this work will mainly extend and apply the framework to odd powers N (Sec. 2), then 
specially discuss the picture resulting for the Airy function when N = 1 (Sec. 3). 

1. Summary of previous results. 

(See [1] for a more detailed review.) 

Definitions and notations: the Schrodinger operator, with all coefficients straight- 
forwardly scaled out, reads as [2-6] 

H d = — 7—x + q N , for N an even integer. (1) 
aq 2 

(Obvious dependences upon the parameter N will be left implied.) 

Over L 2 (R) , H is a strictly positive operator; it has a purely discrete (nondegenerate) 
spectrum {\k}k=o,i,2,..., an d Xk f +oo according to the asymptotic Bohr-Sommerfeld law 

b \£ = 27v(k + 1/2) + o(l), /c^ooinN, (2) 

which uses the classical action 6o^ M (at the energy A, for the Hamiltonian p 2 + q N ), with 

^defiV + 2_ the growth order, and b = -^ r (^)/ r (J + jf)- ( 3 ) 



t Corrections (July 2003) indicated by footnotes. 
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Other dynamical constants describe the phases a,u> such that H\l 2 (u>r) is unitarily equiv- 
alent to aH\L2(^y. this happens for uj = a -1 / 2 and a in the cyclic group {e liip \ £ = 
0, 1,...,L-1}, with 

An N 
<p = — — - : the symmetry angle, L = — + 1 : the symmetry order. (4) 

Since commutes with the parity operator P, it splits as d = H(l ± P); the 
eigenf unctions of positive and negative parity carry the even and odd labels k respectively. 
Spectral functions, defined next, may also be split accordingly. 

The spectral zeta functions: [4,5] 

S(*) = £V (and Z±{s)^ £ A" s ) (Re«>^); (5) 

z- , even 

fc fc odd 

they extend meromorphically to all of C, with only simple poles located at (or among) 
s = (1 — 2j)fx, j = 0, 1, 2, . . .. For Res < \x the series (5) diverge and the definitions must 
be regularized, according to eq.(2); e.g., in each parity sector, 

E V = { E V " ^r^T + k?) if > 

V ^ +ool fc^r 4tt(-s + /U ) 2 k J p ' ( 6 ) 

with both fc, if even for Z + , odd for Z~ . 
In particular, ^ ± (0) = ±1/4 (using eq.(2)), and 

Z ±/ (0) = ^mj- £ lo S^ + ^(logA* - -) - ^logA K } for k,K , (7) 



which in turn generates the (zeta-regularized) determinants, as detH ± = exp — Z ±f (0) 



The spectral (or functional) determinants: [5,7] D ± (X) = f det(H ± + A): 
these are entire functions, more explicitly given by 

D ± (X) = exp[-Z ± '(0)] Yl i 1 + A / A fc) ( for A» < ^ > 2 ) (8) 

, even 
fc odd 

OO r7±/ \ 

" Z± '(°) " E —(-^) n (for a* ^ 1, i.e., N ^ 2, and |A| < A ). (9) 



n 

71=1 



The harmonic oscillator (AT = 2) has a special status in the family: it is the solvable 
case but also the "confluent" case, i.e., the growth order becomes integer (fx = 1), and 
this invalidates several formulae in their generic form given here; in particular, Z ± (l) are 
infinite and eqs.(8-9) diverge, a valid substitute specification being 

D ± (X) = 2 ±1 / 2 v / 2^ 2" A / 2 /r( 2T ^ + A ) for N = 2. (10) 



Hence N = 2, potentially the only elementary example for the formalism under review, is 
a pathological case instead. 

The determinants admit semiclassical asymptotic expansions for | A| — > oo (| arg A| < n — 6). 
These have quite stringent forms, actually simpler for the full determinant D = D + D~ 

and the "skew" ratio D p = D+/D~ : [4] 



logD(A) ~ ^Ta^ 1 " 2 ^ (N^2) (witha = (2sin7r^)- 1 6 ), 

3=0 

oo 

lo gJ D P (A) ~ - l -\ogX + J2 d r X ~ (1+N/2)r - 



r=l 

The leading behaviors imply that D and D ± , as entire functions, are of order ji (the higher 
coefficients a,j, d r are computable term by term, too). 

Exact functional relations: [8,1] the determinants satisfy a basic functional rela- 
tion, obtained by an exact WKB calculation (and corresponding to the reflection formula 
for r(z) when N = 2): 

e^/ 4 D + (X)D-(e^X)-e-^/ 4 D+(e^X)D-(X) = 2i, (N ^ 2) 

i.e., e^/ 4 D+D^ - e"^/ 4 D+D^ = 2i, ' ' 

upon introducing the generic shorthand notation Di(-) d = D(e liip •) (cp being the symmetry 
angle, and t an integer mod L). Within identities like (12), all such subscripts are globally 
shiftable mod L, since A can be freely rotated by e lv throughout. 

Eq.(12) was found [1] as an equivalent form of a multiplicative coboundary formula 
linking the full and skew determinants, [5] 

6_W2 ^P(I) A) = ex Ph 2i$ ( A )]' fOT = arcsin [{D{^X)D{X)Y 1, \ (13) 

D p 

i.e., =expi(-2^ + ^/2), t = 0, 1, . . . , L - 1 (mod L) {N ^ 2). (14) 

This in turn entails a consistency condition upon D, embodied in this cocycle relation (of 
length L), 

L-l 

^$ £ (A)=Lv?/4 =tt/2 for even N (N ^ 2). (15) 

The full determinant D thereby inherits an autonomous functional equation, circularly 
symmetric of order L (and convertible to a polynomial form). 

Spectral sum rules: [5] algebraic identities arise from expanding either functional 
relations (12) or (13) to all powers A n around A = 0, with the help of eq.(9); by eq.(5), the 
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output is an infinite sequence of sum rules connecting various spectral moments: 

[Z'(0)=] Z +/ (0) + Z-'(0) = logsin| (ra=0)(16.0) 
sin | Z+(l) -sin^ Z~(l) =0 (N ^ 2) (n = 1) (16.1) 

sin^ Z+(2)-sin^ Z~(2) = sin £ (2 cos f [Z+(l) - Z"(l)]) 2 (n = 2) (16.2) 

sin(n — i)^ Z + (n) — sin(n + |)f Z~(n) = P n {Z + (m), Z~(m)}i< m<n (rank n) (16. n) 

where P n is a polynomial, homogeneous of degree n under the ruling that each Z ± (m) is 
of degree m. The relative weights in the left-hand-sides recur with period L in n, yielding 
Z(n) whenever n is a multiple of L. For such n = mod L, Z(n) can be further reduced 
to a polynomial in the {Z(m)}i< Tn<n only: this follows by directly expanding the closed 
functional equation for D itself, eq.(15), in place of eq.(12) (an option which however 
misses all the other identities of ranks n ^ mod L). By contrast, no Z + (n) or Z~(n) 
ever comes alone on a left-hand side, which precludes similarly closed functional equations 
for D + or D~ separately. 

Exact quantization conditions: [6,1] exact quantization formulae a la Bohr- 
Sommerfeld for the spectrum ultimately derive from eq.(12), as 

2E + (A fc ) = fc+i + £, fc = 0,2,4,... def iV-2 

l k with K = aTT9 ' (17) 

2E_(A fc ) = fe+---, fc=l,3,5,... 

where, excepting the case N = 2 (singular but solvable otherwise, to S±(A) = A/4), 

E ± (A) d = 7r- 1 ArgD ± (-e- i ^A) t (A > 0), (18') 

ArgD ± (-e- i ^A)= ^ fc/ (A), t </>fc'(A) = f arg(A fc / - e -1 ^ A); (18") 

fc' ev fi 1 
K odd 

i.e., ArgiP 1 * 1 denotes the determination of argD ± which is continuous over the half-line 
[0, — e~ lip oo) t and vanishing at A = 0, and this is the sum of the angles 0fc'(A) (e [0,7r)) 
subtended by the vector (0, a) at the points e 1(p \ k >. 

The system (17-18) appears to specify the exact spectrum self consistently within each 
parity sector. It amounts to a fixed-point condition for the (nonlinear) mapping defined 
by the application of eqs.(18) followed by (17), upon sequences of trial levels A^° that are 
asymptotically correct ([6,1]: to be more precise, we require their compliance with eq.(2)). 

This fixed-point mechanism is moreover easy to implement numerically within some 
preset tolerance e, which has to be nonzero to render the computational scheme finite- 
dimensional (by inducing a high-/c cutoff): upon such truncations, then, the straightfor- 
ward loop iterations of these maps appear to converge towards the correct spectra, consis- 
tently within 0(e), and at geometric rates [6]. Such a constructive enactment of eqs. (17-18) 



t -e- 



e 1Lp and not e 1Lp as in the printed version. 
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works if, essentially, some form of contractivity holds for the underlying map. What we 
now have in this direction is no rigorous statement yet, but the numerical observation of 
such a behavior in every tried case, plus partial analytical evidence: e.g., the mapping 
can be simplified by treating the spectrum as continuous [1], in which approximation its 
contraction factor is just k — and \k\ < 1 by eq.(17). 

2. Extension to odd degrees and applications. 

The above framework will now be strengthened in two ways, both involving a violation 
of analyticity through the invocation of odd degrees N, in which case the confining potential 
becomes Ig]^. 

From even to odd N: all previous results actually follow from solving a connection 
problem for the Schrodinger equation (H + \)ip(q) = along a half-line, e.g., from q = to 
+oo. The required complex WKB calculations were performed exactly [5] thanks mostly 
to the analyticity of the potential (whose homogeneity also eased the operations but at a 
more concrete level). 

Consequently, these procedures apply to odd polynomials q N as well, with the under- 
standing that all results refer to the spectra specified by the same boundary conditions on 
the half-line [0, +oo) as in the case of even N: specifically, at q = 0, Neumann for positive 
parity, resp. Dirichlet for negative parity. But these in turn precisely select the positive, 
resp. negative parity spectra of the potential \q\ N also when N is odd. 

All previous formulae extend to odd N under this interpretation, with just one forced 
explicit change: the order L of the symmetry group, being an integer, cannot retain the 
value N/2 + 1 from eq.(4), but becomes twice that instead: 

L = N + 2 for odd N, (19) 

while (p stays unchanged] consequently, the cocycle relation (15) becomes 

L-l 

®e( x ) = L<p/4 = 7T for odd N. (20) 



• Example: N = 3. We may test the extended formalism upon the confining homogeneous 
cubic oscillator (the operator H = — d 2 /dg 2 + \q\ 3 ), for which 

22/3 /o 

^ = 5/6, </? = 47t/5, L = 5; b = r T(l/3) 3 = a ; k = 1/5. (21) 

Among the several results which can be validated, we describe the exact quantization 
formulae (17-18) as they probe the functional relation (12) the most sharply. So, we 
iterated eqs. (18-17) (under a tolerance value e = 10~ 9 ) upon the input data consisting of 
each parity subsequence of the semiclassical spectrum (given by bo 5//6 = 27v(k + 
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|)), and we readily observed convergence to these lowest eigenvalues: 

k even k odd 

1.0229479 1 3.4505627 

2 6.3702932 3 9.5220764 

4 12.870297 5 16.369373 [ ' 

6 20.000879 7 23.745471 

8 27.592421 9 31.530790 

The convergence appears to be governed by contraction factors +0.233 for positive 
parity (k, k' even) and +0.189 for negative parity (k, k' odd) (cf. the continuous-spectrum 
prediction = k = 1/5). 

Separately, we diagonalized the matrix of H in the orthonormal eigenfunction basis 
{ip n } of the harmonic oscillator — d 2 /dg 2 + q 2 , using 



niN ,, . / g+gV! n"T / l^V r(j(l + Jy + n' + n") - (m' + m")) 

(Vv||<z| Wn»)-\) - 2^ m/Vm"\ (n'-2m')\ {n"-2m")\ 

0< m '<n'/2 V ^ V ' 

0<m"<n"/2 

(23) 

for n' = n" mod 2, zero otherwise; in each parity sector we applied several truncations 
(in the size range 20-40), and retained the eigenvalue figures as far as they were fully 
stable. These results then showed complete agreement with eq.(22). (Here, the brute-force 
calculation might well be the less reliable one: eq.(23) generates numerous large entries, 
through cancellations between even larger terms.) 

The exact method is thus validated to an accuracy of 10 _6 -10 -7 , whereas its input 
(we used the semiclassical spectrum) was off by as much as 10 _1 (for the ground state: 
a[, 0) « 0.920791). 

• General iV: an alternative exact approach exists for odd and even degrees alike: it takes 
the fully analytic potential q N (nonconfining for N odd) and privileges a different spectral 
function, namely a Stokes multiplier — normalized here as C(A) = e -1 '' 3 / 4 c(A) relative to 
refs. [9,2]. The present analysis connects to that approach and recovers its results, through 
the identity 

C = (2i)" 1 (e i ^ 2 D+Dz - e-^/ 2 D+D^) (for N ^ 2). (24) 

For instance, when iV = 3, C(A) obeys a closed functional equation, C2C3 — Co = 1 
(the subscripts being shiftable mod 5; cf. [9]; [2] chap. 5, sec. 27; [10]). Here, this equation 
is not just verified by the expression (24) (using eq.(12)), but even solved thereby, to the 
extent that can be constructed by using the fixed-point eqs. (17-18) to get their zeros 
(cf. eq.(22)) and then forming the Hadamard products (8). (By contrast, we do not know 
how to specify any of the other entire functions, like D or C, as directly). 

(Note: remarkably, the same N = 3 functional-equation structure C2C3 — Co = 1 also 
appears in integrable 2D field-theory models that involve dilogarithms [11], and in the I-st 
Painleve function [12].) 
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Duality: within homogeneous problems, it is worthwhile to ask which analytic po- 
tentials q N , q N might share the same rotation symmetry group, since the latter imprints 
the structure of spectral functions in a basic way. The answer is a duality relation, 

(p + <p' = 2tv 

11 / / ( 25 ) 

implying - + — = 2, L = L , k = and NN = 4. 

One solution is iV = N' = 2: thus the harmonic oscillator is selfdual (but also singular). 
However, there exists one (and only one) other solution: the pair N = 4 (homogeneous 
quartic oscillator [5]) and N' = 1 (Airy equation [13]). 

Hence, the Airy function (as the relevant solution of the linear potential) turns out 
to be conjugated with the spectral determinants of q 4 ! This form of duality may be weak 
(not implying exact links between the solutions of both problems) but still liable to reflect 
strong structural resemblances (now that the framework accommodates even and odd A" 
on the same footing) . 

3. Airy function vs quartic oscillator. 

We will now basically scan properties of the Airy function (N = 1) [13] as they appear 
in the global setting of the \q\ N problem, and especially analyze its duality with the quartic 

case (AT = 4) [4-5]. 

Already, (unlike A" = 2,) both N = 1 and N = 4 share the status of regular values for 
the whole formalism as they give non-integer growth orders fx. Thus, the various constants 
corresponding to N = 4 are 

// = 3/4, v? = 2tt/3, L = 3- b = T(l/4) 2 , a = b /V2; k = 1/3. (26) 

Then (almost) all previous formulae also hold in their raw (generic) form for A^ = 1, upon 
simply resetting all the constants: 

H = 3/2, cp = 4tt/3, L = 3; b Q = 8/3, a = -4/3; « = -1/3. (27) 

Special features of the Airy case: when N = 1, the eigenvalues are those of 
the potential \q\, namely (up to sign) the zeros of the Airy function for negative parity 
and those of its derivative for positive parity (the connection with their usual notation 
[13] being a s = — A2 S -i, a' s = — A2 S -2, s = 1,2,...); the Bohr-Sommerfeld rule (2) 
(plus its corrections, omitted here) reproduces their known asymptotic behavior. The 
corresponding spectral determinants are given by 

D~(X) = 2v^FAi(A), D+(\) = -2 v / ^Ai'(A). (28) 

(Normalization, the only non-obvious issue, is determined by adjusting the asymptotic 
forms (11), which are expressly without constant terms, against the known asymptotic 
expansions of Ai and Ai' [13].) 
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Thus, the extension of the formalism to odd N reveals its only elementary regular 
example, given by the Airy functions (notwithstanding that their underlying potential \q\ 
is the most singular one at 0). 

Still, one amendment is specially mandated when N = 1, namely a (standard) regu- 
larization [7] for those series and products which turn divergent as \x > 1: thus, Z ± (l) are 
now computed from eq.(6) instead of (5), and eq.(8) for D ± accordingly gets replaced by 

D ± {\) = e-^'W+z^Dx JJ (i + A/A fc ) e" A /^ . (29) 

, even 
k odd 

The basic determinantal identities: we first compare the functional relations (12) 
and their consequences for the two conjugate cases, 



e™/ 6 D+(\)D-Q\)-e-™/ 6 D + (j\)D-(\) = 2i for N = 4 {<p = 2tt/3) ^ ^ 

e^/ 3 D + (X)D-(fX)-e-^/ 3 D + (fX)D-(X) = 2i for N = 1 (tp = 4tt/3) J 

(30) 

The latter, by eq.(28), is simply the classic Wronskian relation between Ai(-) and Ai(j 2 -), 
which both solve the same Airy equation: 

W[Ai(.),Ai(j 2 -)] = (27r)- 1 e i -/ 6 . (31) 

Here the Airy function Ai(z) acts in two ways, as solution and as spectral determinant of 
the Airy equation; this confusion of roles stems from the property (specific to N = 1) that 
the operator H + X solely involves the combination variable z = q + X. 

For N > 1, by contrast, the determinants are not known to solve any second-order 
differential equation; while this remains to be confirmed in full generality, the fact is already 
certain for = 2 (due to the presence of F(z): cf. eq.(10), and [2], chap. 5, sec. 27). Still, 
the N = 1 and N = 4 functional relations are highly similar (duality!); but their different 
phase prefactors will suffice to create an essential distinction (as eq.(33) will show). 

If the discrete symmetry rotations are performed upon either of eqs.(30), a closed 
system of L = 3 equations follows: 

e i<^/4 D + D - _ e -i<p/4 D + D - _ 2i 

e W4 D + D - _ e -i<p/4 D + D - _ 2i (3 2 ) 

e i¥>/4 D + D - _ Q -i<p/4 D + D - = 21 

With the DJ as unknowns, this is a linear system, whose 3x3 determinant has the 
value 

A = 2i sin3</V4 D+D+D+. (33) 

• N = 4: if = 2n/3, thus A ^ 0; then the odd spectral determinant D~(= Dq) can be 
solved as a functional of the even one (and vice- versa), in rational terms: 
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(the subscripts being shiftable mod 3). 

• N = 1: (f = 47r/3, thus A = 0, a sharply different situation. Now the system cannot 
be solved for the Dj , and the elimination process yields a (linear) relation among these 
instead, (and likewise for -D + ), 

£-+^-+^£,- = 0, i. e ., Ai(-)+jAi(j-)+j 2 Ai(j 2 -) = 0, 
J D++j J D++j 2 J D+ = 0, i.e, [Ai(.)+jAi(j.)+j 2 Ai(j 2 .)]' = 0; 

hence both identities simply reflect the classic 3-solution dependence relation for the Airy 
equation (but now in partnership with the nonlinear eqs.(34)). 

The cocycle functional equations: their comparison will already uncover a less 
obvious (new?) identity for the Airy functions. 

• N = A: in the quartic case, the cocycle relation (15) (of length 3) implies sin $2 = 
cos($o+$i); expanding and squaring the cos &i away in succession yields 2 sin$ sin $1 sin <E> 2 + 
sin 2 $ + sin 2 $1 + sin 2 $ 2 = 1 and then, by the definition of $ in eq.(13), [5] 

D D 1 D 2 = D + D 1 + D 2 + 2. (36) 

• N = 1: the search for a non-trivial Airy counterpart of eq.(36) must also be directed at 
the full spectral determinant, which now reads as 

D = -4tt Ai Ai' = -27r(Ai 2 )' ; (37) 

the relevant cocycle relation has again length L = 3, but this time it is eq.(20); it now 
implies cos $2 = — cos($o + $1) an d then, proceeding just as before, 

Dl + D\ + Dj - 2(D t D 2 + D 2 D + D D 1 ) + 4 = 
D Q = -27r(Ai 2 )'(-), D x = -27r(Ai 2 )'(j 2 ■), D 2 = -27r(Ai 2 )'(j ■)• ^ 

Thus, (Ai 2 )' exhibits a functional equation which is nonlinear, inhomogeneous, and with 
ternary symmetry like eq.(36) (but rederivable by elementary means, a posteriori) . 

The Stokes multipliers: here, their expression (24) plus the functional relations 
(12), (36) lead to diverging behaviors in the two cases. 

• N = 4: the Stokes multiplier C(A) and the full determinant D(X) get related both ways, 

D = CiCo-l Co = (D D 2 -1) V 2 (39) 

(the subscripts being shiftable mod 3); hence C stands equivalent to D as a spectral 
function, and indeed it has a functional equation very close to eq.(36), [10] 

C0C1C2 =C + Ci + C 2 . (40) 

• N = 1: C(A) degenerates to a constant (= 1) [2] hence it has lost all spectral information. 
(More generally, D(X) might be unrelated to C(A) for odd N.) 
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Therefore, duality is devoid of content under this specific angle. 

The spectral sum rules: 

• N = 4: for the quartic case, the simple substitution <p = 27r/3 in eqs.(16) yields [5]: 

Z'(0) = -log2 (41.0) 
V(l)-Z-(1) = (41.1) 

Z+(2) - V(2) = ?(Z+(1) - Z~{l)f (41.2) 

Z(3) = ^Z(l) 3 - \z{l)Z{2) etc. (41.3) 

The left-hand-side coefficients then recur with period L = 3 in n; all Z(n) with n = 3p 
reduce to polynomials in the {Z(m)}i< m<n only, deducible directly from the functional 
equation (36) for D alone. 

• N = 1: by setting ip = 4n/3 in eqs.(16), similar sum rules are produced for the moments 
of the zeros of Ai(A) for odd parity, resp. Ai'(A) for even parity (i.e., Z~(n) = J2 s (—a s )~ n , 

Z + {n)=Y, s (-<)- n )i 

Z'(0) = - log(2/>/3) [= Z+'(0) + Z-'(0)] (42.0) 
Z+(1) = (42.1) 
Z"(2) = Z-(1) 2 (42.2) 

Z(3) = Z-(l) 3 -|z-(l)Z+(2) = |z(l) 3 -|z(l)Z(2) etc. (42.3) 



All properties of the quartic case persist (L = 3 again) but in addition here, the closure 
property of identities within zeta values of the same parity label holds not only for Z(3p) 
but also for Z + (3p + 1) and Z~(3p + 2), thanks to the autonomous functional equations 
(35) for D + and D~ respectively. 

Now, moreover, the Airy function is also a special function, whose Taylor series is 
accessible by other means [13], and this amounts to a full knowledge of the expansion (9). 
The "spectral" sum rules above thereby get completed by the "special" sum rules below, 

in which p d = - Ai'(0)/ Ai(0) = 3 5 / 6 (27r)" 1 T(2/3) 2 w 0.729011133: 

Z+'(0) - Z-'(0) = - log p (by definition) (43.0) 

Z~(l) = -p (43.1) 

Z+(2) = l/p (43.2) 

Z~(3) = -p 3 + 1/2 (<^>) Z + (3) = 1 (etc.) (43.3) 

This last identity is exceptional for being so simple (the inverse cubes of the zeros of 
Ai'(-z) sum up to unity), and rational: all higher Z ± (n) can be thus expressed when 
N = 1, but as nontrivial (rational) functions of p, hence are never rational-valued again. 
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Eqs. (42-43) remind of known sum rules for integer powers of zeros of the Bessel 
functions J v ([14], and refs. therein) but do not strictly correspond to them, since Ai(z) is 
a Bessel function of a non-integer power of z (and we know of no earlier explicit mention 
of any sum rule for Airy zeros). Both sets of identities however have the same abstract 
basis: they express the consistency between two full specifications of the given function, a 
Taylor series like eq.(9) and a Hadamard product like eqs. (8), (29). 

• For general A > 2, by contrast, the N = 1 explicit values have extensions of the same 
form only for Z ±7 (0) and Z ± (l) ([4], [5] App. C-D: these refs. stated the results for 
N = 2M only [4], but this limitation is unwarranted): 



(0) = log [(N + 2) ^v^/ r(£)] ( 



4n 

' = N + 2 



*♦(!)- z~ m - <V5r> [^] * - \ r (£) r (|) r(g)/r(§ + £ 



(44) 



(the values for other parity labels follow through the identities (16.0-1)). 

The exact quantization of eigenvalues: the exact quantization mechanism (17- 
18) is mainly governed by the symmetry angle ip, and it degenerates, with quantization 
becoming fully explicit, for ip = n (harmonic case). With respect to this "critical" system, 
the quartic oscillator spectrum (with ip = 2-7r/3, k = +1/3) and the Airy zeros (with 
ip = 47r/3, k = —1/3) assume exactly mirror-symmetric positions; in particular, a key 
ingredient of the exact formula, the kernel function arg(A' — e~ lLf> A) (associated with the 
linear "flux operator" [6,1]), takes opposite values in the two cases. 

However, this symmetry is lost at some stage of the procedure (and we see no manifest 
link between the two resulting spectra either). Thus, the system (17-18) is consistent and 
stable under iteration only when applied to spectra having the right growth order /j,, which 
differs in the two cases (3/4 for N = 4, vs. 3/2 for N = 1). These consequences follow: 

• N = 4: since \i < 1, eqs. (17-18) hold in their original form (fig.l, upper part). When 
applied to (definite-parity) input sequences A^ only subject to eq.(2) (i.e., bo A^ 3,/4 ~ 
2n(k + 2), with bo from eq.(26)), the iteration scheme exhibits contractive convergence to 
the exact eigenvalues, with ratios ~ +0.392 for positive parity (k, k' even) and +0.333 for 
negative parity (k, k' odd) [6]. 

• N = 1: now \x > 1 hence the series in (18") diverges like (Asiny?) times the series (5) for 
Z ± (l); its proper regularization, as dictated by eq.(6) for s = 1 and eq.(29), is then t 



Arg J D ± (-e"^A) = lim { V fc ,(A) - Asin^ - A^ /2 | for k', K —fa ■ (45) 



Here the summands </>fc/(A), still defined as in eq.(18") but with <p = 4tc/3, are all negative 
and decreasing (fig.l, lower part); it is now the counterterm = +A(v / 3/7t)aJ/ 2 which out- 
weighs the whole sum to produce a positive and increasing left-hand side, as required for 
eq.(17). 



t See footnote at eqs. (18). 
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Thereupon, the numerical scheme behaves for N = 1 as for the quartic case. When 
applied to each parity subsequence of the semiclassical spectrum (given by §a£ S ^ 2 = 
27r(/c + i)), the iteration of eqs.(17),(18'),(45) exhibits convergence to the Airy zeros, with 
(roughly estimated) contraction factors —0.37 for positive parity (zeros of Ai'; k, k' even) 
and —0.25 for negative parity (zeros of Ai; k, k' odd). (However, whereas for A > 2 
this scheme is comparatively quite efficient, for N = 1 it should not beat special-purpose 
algorithms for obtaining the Airy zeros; also in this case, the regularization (45) proceeds 
through large cancellations which strongly degrade the final accuracy.) 

In conclusion, exact fixed-point quantization works for the spectra of Airy zeros too. 
Then, moreover, the conjugation symmetry shows this spectrum quantization to be exactly 
as nontrivial as that of the quartic oscillator, even though the Airy equation itself is more 
elementary than its quartic — or even its harmonic (and trivially quantized) — analog. 
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Fig. 1. Geometrical depiction of the exact-quantization summands 0fc'(A) (at A = 15.), in 
the odd spectral sector. Upper half-plane: for the quartic oscillator, using (p = 2n/3 in 
eq.(18"); lower half-plane: for the zeros of Ai(— A), using (p = in eq.(45). 
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